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Recently Carr and Wu (2004, 2005) and also Huang and Wu (2004) show that most stochastic 
processes used in traditional option pricing models can be cast as special cases of time-changed Levy 
processes. In particular these are models which can be tailored to exhibit correlated jumps in both the log 
price of assets and the instantaneous volatility. Naturally similar to a recent work of Barndorff-Nielsen 
and Shephard (2001a, b), such models may be used in a likelihood based framework. These likelihoods 
are based on the unobserved integrated volatility, rather than the instantaneous volatility. James (2005) 
establishes general results for the likelihood and estimation of a large class of such models which include 
possible leverage effects. In this note we show that exact expressions for likelihood models based on 
generalizations of Huang and Wu (2004) and Carr and Wu (2005), follow essentially from the arguments 
in Theorem 5.1 in James (2005) with some slight modification. We show that that an explicit likelihood 
for any of these types of models only requires knowledge of the characteristic functional of a suitably 
defined linear functional. This serves to formally verify a claim made by James (2005). 


1 Introduction 

Recently Carr and Wu (2004, 2005) and also Huang and Wu (2004) show that most stochastic 
processes used in traditional option pricing models can be cast as special cases of time-changed Levy 
processes. This includes for instance the models of Duffie, Pan, and Singleton (2000) and leverage 
effects model of Barndorff-Nielsen and Shephard (2001a, b). In particular these are models which 
can be tailored to exhibit correlated jumps in both the log price of assets and the instantaneous 
volatility. Naturally similar to a recent work of Barndorff-Nielsen and Shephard (2001a, b), such 
models may be used in a likelihood based framework. James (2005) establishes general results for 
the likelihood and estimation of a large class of such models, based on quite general linear functionals 
of Poisson random measures, which include possible leverage effects. In this note we show that exact 
expressions for likelihood models based on generalizations of Carr and Wu (2005) and Huang and 
Wu (2004), follow essentially from the arguments in Theorem 5.1 in James (2005) with some slight 
modification. This serves to formally verify a claim made by James (2005). We shall be rather brief 
in our exposition and refer the reader to the above mentioned works for further references and 
motivation. Huang and Wu (2004) and Carr and Wu (2005) proposed a model for the log price of 
assets which can be written as 

(1) x*{t) = (r - q)t + Jiird)) + J 2 {l{t)) + PtU) + aj(t) + (TlTi(r(f)) -k 

where (Ji, J 2 ) are independent pure jump Levy processes,(Wi, W 2 ) are independent standard Brow¬ 
nian motion independent of (Ji, J 2 )- Furthermore (r, 7 ) are non-negative random time changes, inde¬ 
pendent of the above processes. The independence property can be relaxed. See Carr and Wu (2004). 
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Integrated Volatility 


An example of r and 7 are the integrated Ornstein-Uhlenbeck models of Barndorff-Nielsen and Shep¬ 
hard (2001a, b) which are used to model the integrated stochastic volatility. For notational conve¬ 
nience we shall hereafter set cr = 1. Assuming conditional independence across intervals [(i — 1 )A, zA] 
for z = 1 ,..., n and A > 0, define := r((zA)) — r(((z — 1 )A)) and 7 ^ := 7 ((zA)) — 7 ((z — 1 )A)), 
additionally define 

Ji,i = Ji(r((zA))) - Ji(t((z - 1 )A)) 

and 

^2 = J2(7((zA)))-J2(7((z-1)A)). 

It follows from o that one may define a likelihood model for the aggregate returns Xi = x*{iA) — 
a;*((z— 1 )A) which is conditionally Normal given (Ji, J 2 , t, 7 ). Assuming that these quantities depend 
on some unknown Euclidean parameter 9, this is in essence the type of framework addressed in 
James (2005) except for the appearance of the time-changed (Ji, J 2 ). Before we proceed to derive 
an exact form of this present model, we will further generalize o by modeling the t and 7 as linear 
functionals of a quite arbitrary Poisson random measure, say N, on a Polish space Y. In this way 
the t in the model o may be replaced by a more abstract notion of time. Note that we can always 
choose A on a big enough space so that r and 7 are either independent or dependent. 

Remark 1. Note that we take the quite general Poisson framework for explicit flexible con¬ 
creteness. This includes both discrete and continuous quantities, such as shot-noise processes on 
abstract spaces. As we shall see, on a more abstract level we simply can specify t and 7 so that one 
knows the explicit characteristic functional of linear combinations of such processes. 


2 Exact expression for the marginal likelihood 

First as in James (2005) let N denote a Poisson random measure on some Polish space 'V with mean 
intensity, 

E[N {dx)\i'] = v{dx). 

We denote the Poisson law of N with intensity v as F{dN\i'). The Laplace functional for N is defined 
as 

]g|'e-'^(/)|jy] = f e“'^^-^^P(dA|zz) = 

Jm 

where for any positive /, N{f) = f{x)N{dx) and A(/) = /^(l — e~^^^'>)v{dx). M denotes the 
space of boundedly finite measures on 'V. We suppose that r = N{hi), for z = l,...,7 z where 
hi,... ,hn are positive measureable functions on 'f'. Similarly we suppose that for each z 7 i = N{gi) 
where gi are positive measureable functions on 'V. With this specification it follows from o that the 
generalized notion of aggregate returns Xi\Ti,yi, Jty, fi,r,q are conditionally independent 

Normal random variables expressible as 


( 2 ) 


Xi = {r - q)A + -I- Jj ,2 + Pin + 7i) + {a - /3)7i -I- 


where are independent standard Normal random variables. Now for fi{x) := hi(x) + gi{x) on Y, 
set 

T* ■■= Ti+y, = N{fi) 

and set n = {r — q) and a — (3 = p. Then the conditional Normal density of each Xi, say 4>{Xi\p.A -|- 
Jt,i + Ji ,2 + Pli + (3 t* ,Ti), can be written as 


g(Ai-Ji,i-Ji,2-P7i)/3 


1 1 


[e-(A-Ji.l-Ji.2-P7i)V(2T;)jg-T*/3V2 


( 3 ) 
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where Ai = Xi — fj,A. It is not difficult to see that if Ji^i and Ji^2 were removed, the subsequent 
likelihood model is a special case of the models handled by James (2005, Theorem 5.1). James (2005) 
remarks that the models such as © pose no additional difficulties as long as one knows the char¬ 
acteristic functional of say Ji i and 2 and their sums. However Ji and J 2 are pure jump Levy 
processes and hence by way of the known Levy -Khinchine formula and independence relative to r 
and 7 , this point is a trivial matter. In particular for any real or complex number, u write for j=l, 2, 

where the explicit form of ipj (w) is given by the Levy-Khinchine formula which can be found in Carr 
and Wu (2005), but is otherwise similar to the function A. It then follows that for each i 

and E[e-“'^*'"] = 

Note also conditional on N, for possibly complex valued numbers (wi,..., w„), 

n n n n 

(4) E[J| 1 ^] ^ JJ- J,, 21 ^] ^ JJ- g-7.V'2(^d, 

i—1 i—1 i—1 i—1 

How we shall proceed is to first evaluate everything conditional on N. Our results will then boil 
down to expectation of exponential sum of terms of the form, 

(5) Ti[ipi{l3 + + {0^ + y‘i)/2] + 7 i['(/' 2(/3 + ^Vi) + p{P + £.yi) + + y^)/^] 

where ^ is the imaginary number. This is similar to the case of James (2005). Now for real valued 
numbers (yi,... ,y„), set 


^n{x) = + ^Vi) + {0^ +yl)l‘Ahi{x) 


and 

n 

"^"(a:) = Xl[^ 2(/3 + (,y^) + p{P + ^y^) + {0^ + yf)l2]g,{x) 

Note that the sum over the terms in is equivalent in distribution to N0ln F T„). We now state 
the form of the likelihood. 

Theorem 2.1 Suppose that N is a Poisson random measure with intensity v on Y. Furthermore 
suppose that Ti and 7 ^, defined above, are chosen such that A(0„-|-T„) < 00 . Then the joint marginal 
density or likelihood of Xi, ..., Xn\fi, P,9, p, determined by ^ and m is given by, 

„nA0 r ” 

^{X\p, P, 0, p) = n 

Jb" 

Where A = 000^\Ai/n. The result applies for the case where N is not necessarily Poisson, by 
replacing with □ 


Proof. The proof of this result is simply a slight variation of Theorem 5.1 in James (2005). 
For completeness we give many of the same details. Here we use the fact that for each i one has 
the identity deduced from the characteristic function of a Normal distribution, with mean 0 and 
variance l/rf, evaluated at Wi = Ai — Ji^i — Ji ^2 — PJi- That is. 
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Integrated Volatility 


Now the result proceeds by substituting this expression in m and applying Fubini’s theorem . One 
then integrates out the expressions involving (Jip, >^, 2 ) conditionally on N, which results in using 0). 
The rest now is precisely as the proof of Theorem 5.1 in James (2005). That is after rearranging 
terms it remains to calculate the expectation of □ 

Remark 2. Note that these arguments may be used to directly evaluate the density of x*{t) 
given observations Xi,... ,Xn, which might be interesting in an option pricing context. Statistical 
estimation follows along the lines of James (2005). 
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